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A  SCALING  AND  PLOTTING  ROUTINE  FOR  TWO  DIMENSIONAL  DATA 


Introduction 


The  recent  advent  of  more  reliable,  more  accurate  and  faster  plotting 
and  graphics  output  devices  and  recent  attempts  by  computer  manufacturers 
to  standardize  graphics  plotting  software  (Core  System^  graphics 
standard,  HP  Graphics-1000,  etc.)  have  required  new  application  software 
development  of  general  scientific  data-plotting  routines.  Our  specific 
needs  in  chemical  experimentation  require  the  capability  for  plotting  of 
a  wide  variety  of  X  vs.  Y  data  types  and  magnitudes,  with  and  without 
"error  bar"-type  error  limits  on  each  data  point.  In  addition,  it  is 
often  necessary  to  do  linear  least  squares  calculations  on  these  same 
data  and  plot  a  least  squares  regression  line  on  the  same  plot  for 
visual  indication  of  linearity,  scatter  and  goodness-of-fit.  Finally, 
the  capability  to  provide  journal-ready  plots  to  eliminate  the  need  for 
the  user  to  make  decisions  about  scaling  is  highly  desirable. 

The  program  GPLOT  satisfies  these  basic  requirements  and  contains 

such  additional  capabilities  as  multiple  data  sets  on  a  single  graph, 

multi-colored  plots  and  variable  origin  starting  location.  The  main 

routine  does  the  plotting;  two  subroutines  are  used  to  scale  the  data 

and  calculate  a  least  squares  regression.  Ail  three  routines  are  written 
(2) 

in  Fortran  IV  .  This  software  was  developed  to  provide  maximum  flexibility 
with  a  minimum  of  effort  by  the  user.  It  is  designed  to  run  on  a  Hewlett- 
Packard  1000  computer  system  under  an  RTE  IV  operating  system  and  makes 
use  of  the  device-independent  features  of  Hewlett-Packard's  Graphics-1000 
Software  as  well  as  some  H.P.  extensions  to  standard  Fortran  IV. 

This  document  is  designed  to  serve  the  purposes  of  a  users'  guide 
and  an  operations  manual,  and  to  provide  sufficient  documentation  for  program 
maintenance. 

Manuscript  submitted  July  31, 1980. 
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Features  of  Program 


1.  Automatic  scaling. 

2.  Least  squares  line  with  slope,  Y -intercept  and  standard  deviations 
of  each. 

3.  Error  bars. 

4.  Multiple  plots  with  different  symbols  on  same  graph. 

5.  Data  source  can  be  from  disc  file,  cartridge  tape  or  typed  in 
from  the  keyboard  of  a  terminal. 

6.  Output  can  be  on  graphics  terminal  or  plotter. 

7.  Plots  can  be  line,  symbols,  or  symbols  connected  by  line. 

8.  Axis  labeling  and  title  are  entered  from  terminal. 

9.  Origin  can  be  determined  automatically  or  forced  to  start  at 

(0,0). 

10.  Tick  marks  are  labeled. 

11.  Multiple  color  plots  with  automatic  pen  changing. 

12.  A  manual  scaling  option  can  plot  multiple  data  sets,  which  have 
different  maximum  and  minimum  values,  on  a  single  set  of  axes. 

13.  The  routine  can  handle  up  to  100  points  per  data  set. 


Computer  and  System  Configuration 

This  program  could  be  modified  to  run  on  many  different  computers 
with  different  plot  packages;  however,  this  routine  was  written  for  the 
Hewlett-Packard  HP  1000  family  of  computers  using  Graphics-1000  (HP 
92840A  graphics  plotting  software).  The  operating  system  is  RTE  IVB 
with  updated  software  revision  code  2001.  A  graphics  device,  HP  2648A 
graphics  terminal  and/or  HP  9872B  plotter,  is  needed  to  do  the  plotting. 
Older  versions  of  the  RTE  operating  system  and  earlier  plotters  can  be 
used  so  long  as  graphics  1000  limitations  are  satisfied.  EMA  and  spooling 
features  were  not  utilized.  This  routine  requires  a  24K-word  partition 
to  run. 


User  Changes 

This  program  was  tested  on  a  system  probably  configured  differently 
from  that  of  the  user.  It  is  the  responsibility  of  the  user  to  make 
changes  to  logical  unit  assignments  to  implement  this  program  on  his 
system.  Changes  will  need  to  be  made  to  lines  21  through  26,  31,  226 
and  294  of  program  GPLOT.  LUG  is  the  logical  unit  number  of  the  graphics 
device.  ID  is  the  identification  number  assigned  by  the  device  link 
table  to  a  graphics  logical  device.  Line  44  of  program  GPLOT  does  not 
permit  a  logical  unit  number  greater  than  30.  This  limit  was  set  to 
avoid  input  errors  and  may  need  to  be  changed  by  the  user. 

Program  Background 

The  Hewlett-Packard  plot  package  (HP  92840A  graphics  plotting 
software)  does  not  contain  a  scaling  routine.  A  major  reason  for  development 
of  this  routine  was  the  absence  of  useable  existing  scale  routines.  Most 
CalComp-type  routines  provide  very  limited  plotting  capabilities.  The 
HP  CalComp  scale  routine  requires  placement  of  a  tick  mark  at  one  inch 
intervals  on  each  axis;  this  requirement  is  very  restrictive  in  maximizing 
plot  size.  A  more  generally  applicable  routine  was  required,  which  would 
provide  a  minimum  of  unused  plotting  area  for  a  wide  variety  of  data  types. 


Cautions 


The  routine  GPLOT  and  the  Graphics-1000  routines  which  it  calls  are 
not  omniscient.  It  is  possible  to  enter  responses  or  data  which  will 
produce  unpredictable  results.  It  is  also  possible  to  use  portions 
of  the  program  or  to  modify  it  in  such  a  way  as  to  produce  other-than- 
desirable  results. 

The  subroutine  SCALE  transforms  the  array  X  into  exponential  notation, 
storing  the  exponent  in  IEXP.  This  process  changes  the  array  X  by  some 
factor  of  10.  If  the  user  calls  this  subroutine  from  any  other  program 
he  should  be  aware  that  any  value  passed  back  to  the  calling  routine  could 
be  changed.  For  this  reason  X  and  IEXP  must  be  used  together. 

Certain  data  characteristics,  such  as  failing  to  separate  the  data  on 
an  input  line  with  commas  or  spaces,  will  not  be  detected  as  an  error  by 
this  program. 

Some  errors  in  input  data  can  generate  error  conditions  in  Hewlett- 
Packard  library  routines.  As  written,  the  program  allows  these  to  be 
printed  on  the  standard  list  device  (logical  unit  6).  To  avoid  this,  the 
user  must  supply  his  own  error  routine  as  described  in  ER0.E^, 

Input 

Data  can  be  input  to  the  program  from  a  disc  file,  cartridge  or 

paper  tape  or  typed  in  from  the  keyboard  of  a  terminal.  Disc  files  must 

be  type  3  or  4^.  The  program  can  handle  an  array  up  to  100  data  points. 

The  format  is  X,  Y,  DELX,  DELY  using  free  field  input.  The  X  and  Y 

pair  is  the  position  of  the  point  along  the  X  and  Y  axis  respectively. 

The  optional  pair  DELX  and  DELY  is  the  standard  deviation  in  X  and  Y. 

These  values  are  used  for  drawing  error  bars  and  will  be  doubled  and  scaled 

to  plotter  units  to  provide  the  horizontal  or  vertical  separation  of  error  bars. 

—38  38 

All  numeric  input  data  must  be  in  the  range  of  10  to  10  . 


Loading  Procedure  on  HP  1000 

The  loader  must  be  loaded  as  a  large  background  program  (type  4)  and 
will  require  a  size  of  24  pages  to  accommodate  the  Graphics-1000  routines. 
After  compiling  &GPL0T,  &SCALE  and  &LSREG ,  execute  the  loader  interactively 
as  follows: 

RU.LOADR, , , ,LB 

RE , ZGPLOT 

RE.XSCALE 

RE.XLSREG 

RE.XDLTBL 

SE.XGPS 

END 

where  XGPS  is  the  Graphics-1000  library  file  created  when  Graphics-1000 
(3) 

was  loaded.  Alternatively,  the  source  version  available  from  the  author 

contains  a  loader  command  file  which  can  be  used.  The  loaded  program  will 

occupy  approximately  23K  words  of  memory.  It  would  be  possible  to  decrease 

(43 

the  main  program  space  somewhat  by  using  EMA  for  array  space.  However, 
the  bulk  of  this  23K  words  is  required  for  Graphics-1000  routines. 

Program  Source  Availability 

A  program  source  is  available  from  the  author  on  a  user-supplied 
Hewlett-Packard  264X-type  cartridge  tape.  This  tape  contains  5  files. 

The  first  file  is  a  description  of  what  is  on  the  tape.  The  second  file 
is  a  command  file  which  may  be  used  by  the  loader  to  load  the  programs. 

The  third  file  is  the  source  for  program  GPLOT.  The  fourth  file  is  the 
source  for  subroutine  SCALE.  The  fifth  file  is  the  source  for  subroutine 
LSREG. 


Error  Messages 

Program  GPLOT  checks  for  errors  that  could  occur  when  reading  a  disc 
resident  file.  Messages  to  the  user  are  sent  with  the  name  of  the  file 
if  a  problem  is  encountered.  The  program  has  several  built  in  checks 
to  catch  typing  errors  by  the  user. 


Testing  of  Program 

The  author  has  used  a  wide  variety  of  data  types,  multiple  files, 
and  several  Input  sources  to  debug  this  program.  In  addition,  the  program 
has  been  used  extensively  by  six  individuals  with  different  applications 
and  occasionally  by  about  twenty  others. 

Flowcharts 


Figures  1  and  2  show  the  logic  flow  in  the  main  program  GPLOT  and  in 
the  subroutine  SCALE.  The  calculations  performed  by  subroutine  LSREG  appear 
in  the  Formulas  section  of  this  report.  A  complete  program  listing  appears 
in  the  Appendix. 

Examples 

Figures  3,  4,  and  5  are  examples  of  different  types  of  plots  that 
GPLOT  can  produce.  Figure  3  is  a  symbol-only  plot  with  error  bars  and  a  lease 
squares  line.  Figure  4  is  a  line  plot  of  51  data  points.  Figure  5  is  a 
plot  of  3  data  sets  using  different  symbols  for  each  data  set.  These  three 
figures  demonstrate  the  minimization  of  unused  plot  area  resulting  from 
the  algorithm  in  SCALE. 

Execution  of  Program 

The  execution  of  GPLOT  can  be  performed  by  supplying  terminal  and 
graphics  device  Logical  Unit  Numbers  (LU  and  LUG)  via  the  run  string 
parameters.  These  parameters  (globals)  are  retrieved  by  a  call  to  the 
Hewlett-Packard  routine  RMPAR  in  GPLOT.  If  these  parameters  are  not 
supplied,  the  program  retrieves  LU  and  prompts  the  user  for  the  graphics 
device  LUG.  Device  selection  is  to  be  made  between  the  plotter,  for  a 
hard  copy,  or  the  graphics  terminal.  The  data  source  is  from  a  disc 
file,  cartridge  or  paper  tape  or  entered  via  the  keyboard  of  a  terminal. 

If  the  source  is  a  disc  file  the  name  of  the  file  is  requested.  If  the 
source  is  other  than  the  disc  a  logical  unit  number  is  requested.  One 


of  three  types  of  plots  can  be  selected;  a  straight  line  connecting  each 
point,  a  symbol  at  each  point  or  a  symbol  at  each  point  with  connecting 
lines.  Labels  are  entered  from  the  users  terminal,  any  ASCII  character 
is  permitted  (capital  and  lower  case  letters,  numbers  and  symbols).  The 
X  and  Y  axis  labels  cannot  exceed  30  characters  and  the  title  of  the  plot 
cannot  exceed  40  characters  in  length;  characters  beyond  these  limits  are 
ignored.  The  starting  position  of  the  origin  can  be  forced  to  start  at 
X  *  0  and  Y  ■  0  or  the  user  can  let  the  scale  routine  determine  an  origin 
that  will  maximize  the  size  of  the  plot  vs.  the  size  of  the  axes.  Scaling 
is  performed  automatically  by  the  SCALE  subroutine,  however  a  manual  override 
is  provided.  The  manual  scaling  mode  is  used  to  increase  the  limits  between 
the  maximum  and  minimum  values  of  an  axis.  This  feature  is  necessary  when 
plotting  multiple  plots  on  one  graph  when  the  maximum  and  minimum  values  of 
all  of  the  data  secs  are  not  within  one  data  set.  When  plotting  multiple 
data  sets  the  first  set  plotted  must  have  the  smallest  and  largest  values  of 
all  of  the  data  sets.  If  this  condition  cannot  be  met  the  user  must  specify 
manual  scaling  and  enter  minimum  and  maximum  values  of  the  entire  set  of  data. 

A  least  squares  line  can  be  drawn  on  the  plot  with  slope,  y- intercept  and 
respective  standard  deviations  printed  on  the  users  terminal  or  printer.  Error 
bars  can  be  drawn  around  each  point  provided  requirements  in  the  input  section 
of  this  paper  have  been  met.  Error  bars  that  are  small  enough  to  distort 
the  symbol  printed  at  a  data  point  are  suppressed  and  a  message  is  printed 
on  the  users  terminal  of  this  action.  Six  or  less  plots  can  be  made 
using  the  same  set  of  axes  provided  all  data  sets  fall  within  the  limits 
of  the  first  data  set  plotted.  The  user  has  the  option  to  make  pen  color 
changes  when  doing  multiple  plots  on  the  same  axis. 

Subroutines 

The  subroutine  SCALE  uses  a  table  look-up  method,  based  on  the  difference 
between  maximum  and  minimum  values,  to  determine  axis  scaling  and  number 
of  tick  marks  to  be  placed  on  each  axis.  The  SCALE  routine  uses  an  algorithm 
that  shifts  decimal  points  to  Increase  numbers  that  are  less  than  one  and 


decrease  numbers  that  are  greater  than  1000.  This  method  can  handle  a 
difference  of  any  order  of  magnitude  and  the  data  trill  be  scaled  to  cover 
a  minimum  of  50 %  of  each  axis. 

The  subroutine  LSREG  does  a  least  squares  linear  regression  calculation 
including  the  standard  deviation  of  the  slope  and  y  intercept. 

Plotter  Setup 

The  HP  98728  plotter  is  used  when  hardcopy  results  of  GPLOT  are 
required.  The  plot  generated  by  GPLOT  is  designed  to  fit  on  standard 
84  X  11  inch  paper  with  adequate  margins  for  publication.  In  order  to 
center  the  plot  the  paper  must  be  placed  \  inch  from  the  left  side  and  1 
inch  from  the  bottom  lower  left  corner  of  the  plotter  bed.  This  displace¬ 
ment  is  necessary  because  the  HP  9872B  plotter  with  advance  option  OFF 
places  the  lower  left  corner  of  the  plotting  window  at  (520,  380)  instead 
of  (0,0).  For  multicolored  plots  pen  placement  is  as  follows;  pen  1  is 
black,  2  is  red,  3  is  green,  and  4  is  blue.  If  the  plot  Is  to  be  one  color 
that  color  pen  must  be  in  pen  holder  1. 

Plotting  Accuracy 

The  HP  9872B  plotter  is  divided  into  plotter  units  where  one  unit  * 
0.025  mm.  This  is  to  say  that  the  overall  resolution  of  the  plot  is  one 
part  in  N  where  N  is  the  number  of  plotter  units  occupied  in  the  X  or  Y 
direction.  The  X  axis  is  7.7  inches  (7,823  plotter  units)  and  the  Y  axis 
is  5.25  inches  (5,334  plotter  units).  The  plotter  resolution  is  one  part 
in  7,823  in  the  X  direction  and  one  part  in  5,334  in  the  Y  direction.  In 
the  worst  case  of  plot  coverage  vs.  axis  length  (50%)  the  plot  resolution 
would  be  one  half  of  the  plotter  resolution. 


Mnemonic  List 


..  »w-  epifpiiu  i„  iypp.  .f.w 


BARX  -  Height  of  X  error  bar  tick  mark 

BARY  -  Height  of  Y  error  bar  tick  mark 

DELX  -  Experimental  error  in  X  (standard  deviation) 

DELY  -  Experimental  error  in  Y  (standard  deviation) 

DIF  -  Difference  between  max  and  min 
HH  -  Half  height  in  character  cells 
HW  -  Half  width  in  character  cells 
IAX  -  Label  of  X  axis  (30  characters  max) 

IAY  -  Label  of  Y  axis  (30  characters  max) 

IBAR  -  Type  of  error  bars 

ID  -  Identification  number 

IDONE  -  Check  for  termination  or  multiple  plot 

IEXP  -  Exponent  of  base  10  in  X  data  scale 

IFMT  -  Source  of  data 

IHED  -  Title  of  plot  (40  characters  max) 

ILINE  -  Check  for  least  squares  line 
IPEN  -  Pen  number  of  plotter 

ISCAL  -  Set  to  zero  for  automatic  scaling,  1  for  manual 

IZERO  -  Set  to  zero  to  force  origin  to  start  at  (0,0) 

JCHAR  -  Character  to  be  plotted  at  each  data  point 

JEXP  -  Exponent  of  base  10  in  Y  data  scale 

JJ  -  Type  of  plot 

KK  -  Number  of  plots  on  same  axis 

LU  -  Logical  unit 

LUG  -  Logical  unit  of  graphics  device 

LUT  -  Logical  unit  of  cartridge  tape  or  keyboard 

NAME  -  Name  of  data  file 

NOBAR  -  Check  to  see  if  error  bar  was  too  small  to  plot 
NP  -  Number  of  points 

51  -  Standard  deviation  of  slope 

52  -  Standard  deviation  of  Y  intercept 
SLOPE  -  Slope  of  least  squares  line 
SXTIC  -  Interval  between  X  tick  marks 
SYTIC  -  Interval  between  Y  tick  marks 
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X  -  Displacement  along  X  axis 

XBAR1  -  Distance  to  right  of  character  of  X  error  bar 
XBAR2  -  Distance  to  left  of  character  of  X  error  bar 
XEND  -  X  value  at  end  of  least  squares  line 
XMAX  -  Maximum  value  of  X 
XMIN  -  Minimum  value  of  X 

XST  -  X  value  at  start  of  least  squares  line 
XTIC  -  Number  of  tick  marks  on  X  axis 
Y  -  Displacement  along  Y  axis 
YBAR1  -  Distance  above  character  of  Y  error  bar 
YBAR2  -  Distance  below  character  of  Y  error  bar 
YEND  -  Y  value  at  end  of  least  squares  line 
YINT  -  Y  intercept  of  least  squares  line 
YMAX  -  Maximum  value  of  Y 
YMIN  -  Minimum  value  of  Y 

YST  -  Y  value  at  start  of  least  squares  line 
YT1C  -  Number  of  tick  marks  on  Y  axis 
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Formulas 


The  following  formulas  were  used  in  subroutine  LSREG: 


Slope  -  (DIF) (IXIY)  -  (IXIY) 
(DIF) (IX) 2  -  (IX)2 


Y-Intercept 


Standard  Deviation 
of  Slope 


Standard  Deviation 
of  Y  Intercept 


IY  -  (SLOPE) (IX) 
DIF 


A 


/ 

\2 

( (IXIY)  - 

IXIY  ] 

\  -  (I'D 

-  V 

DIF/ 

DIF 

ix2  - 

imi 

DIF 

(DIF 

-  2) (IX2)  - 

(IX)2 

DIF 

2  7Y* 

y\\  (SD  of  SLOPE} 


DIF  ■  the  interval  over  which  the  calculation  is  computed  (FR0M-T0+1) 
IX  *  summation  of  X  values 
IY  »  summation  of  Y  values 


Fig.  2  —  Program  scale 
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INPUT  FOR  FIGURE  4 


391.32234 

397.41956 

403. 31677 

409 . 2i39V 

415 . 11113 

421 .00636 

426.90558 

432.80280 

438.70001 

444.59723 

450 . 49445 

456 . 39166 

462.28888 

468 .18610 

474.08331 

479.98053 

485.87775 

491  .77496 

497.67218 

503.56940 

509.46655 

515. 3o377 

521.26099 

527. 15820 

533.05542 

538.95264 

544 . 34985 

SSQ. 74707 

556 . 64429 

562. 541S0 

568.43872 

574 . 33594 

S80 .23315 

586 . 13037 

S92. 02759 

597.92480 

603.82202 

609.71924 

615.61646 

621 .51367 

627.41039 

633. 30311 

639.20532 

645.10254 

650 . 99976 

656 . 89697 

662.79419 

668.69128 

674.53850 

680 . 48572 

686.33293 
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. 34266162 

.44373513 

. 57228398 

. 69693267 

. 77269948 

.81021631 

.82622504 

.81192720 

.77404330 

.70316505 

.62104368 

.52572405 

.43347001 

. 35731497 

.29903460 

.26176220 

.24086520 

.2314':'  43 

.  2297*4453 

.23414394 

.24093742 

. 24440914 

.25149703 

.26249540 

.27337164 

.28766954 

.31137723 

.33948427 

.36913002 

. 40791887 

. 43468165 

.45912260 

.47635341 

. 48246366 

.  47574240 

. 45777333 

. 43223757 

.39129901 

, 34364962 

. 30037386 

.25235242 

. 19369423 

.  14542344 

.10436270 

.  723451 08E-01 

.  5181473SE-01 

.  287 1 6 073E-0 1 

. I 625320 7E- 01 

. 91 653429E-02 

. 4765 978 IE- 02 

.  4765?78lE-02 


33NVS0OS0V 


INPUT  DATA  hQR  FIGURE  S 


PLOT  i 

5810.32.4 

9140.127.3 

12240 .163.5 
15330,20 9. 1 

19130.256.3 


PL  Of  2 

8450.97.2 

11720 . 130 . 2 
14690 , 167. 1 

19360.199.6 

10810.231.6 


PLOT  3 

9108.138.3 
11500,220.1 

16040.253.3 
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TIME  <MICRO  SEC.  >  X  10** 

Fig.  5  —  Sample  #  R3267 
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Appendix  A 

-  SAMPLE  DIALOGUE  OP  PROGRAM  EXECUTION 


; RU>  GP LOT 

OUT  RUT  ON  GRAPHICS  TERMINAL  TYPE  0 
ON  PLOTTER  TYPE  i 

1 

*#*  SOURCE  OF  DATA  *** 

DISC  FILE  TYPE  0 

KEYBOADRD  OR  TAPE  TYPE  i 

0 

ENTER  NAME  OK  DATA  FILE 
TKOi 

***  TYPE  OF  PLOT  **# 

LINE  PLOT  TYPE  i 

SYMBOLS  CONNECTED  WITH  LINES  TYPE  2 

SYMBOLS  PLOT  TYPE  3 

X 

ENTER  X-LABEL , Y-LABEL  4  TITLE  ON  3  SEPARATE  LINES 
WAVELENGTH  (nn) 

ABSORBANCE 

FIGURE  4  !  PLOT  OF  A  UV  SPECTRUM 
**#  ORIGIN  LOCATION  *** 

TO  FORCE  ORIGIN  TO  START  AT  (0,0)  TYPE  0 
TO  LET  SCALE  DETERMINE  ORIGIN  TYPE  1 

i 

*##  SCALING  #*# 

AUTOMATIC  TYPE  0 

MANUAL  TYPE  1 

0 

LEAST  SQUARES  LINE  ? 

NO  TYPE  0 
YES  TYPE  1 

0 

ERROR  BARS  ? 

NONE  TYPE  0 

X  TYPE  l 

Y  TYPE  2 

X  4  Y  TYPE  3 

0 

TO  EXIT  TYPE  0 

IF  YOU  WANT  ANOTHER  PLOT  ON  SAME  AXIS 

AND  ALL  X  4  Y  VALUES  ARE  WITHIN  THE 

SCALES  OF  THE  FIRST  PLOT  TYPE  I 

0 
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Appendix  B 

r_.1A  ,  PROGRAM  LISTINGS 

i-  T  N4  ,  L 

PROGRAM  GPLOT'  4  ,v9) 

C  GENERALIZED  GRAPH  PLOTTING  ROUTINE  FOR  THE  HP  9372B  PLO'iTER 

C  MAXIMUM  100  PTG .  ALLOWED,  X  AND  Y  MAY  BE  INPUT  FROM  A  DISC 

C  FILE,  CARTRIDGE  TAPE  OR  KEYBOARD  FORMATTED  IN  <X,Y>  PAIRS 

C  FOR  EXPERIMENTAL  ERROR  BARS  FORMAT  IS  ( X , Y , DELTAX , DELTAY ) . 

C 

DIMENSION  1DCBU44)  ,  NAME  (3)  ,  IAX(iS)  ,  IAY(i5) ,  I  HE'D  <20  >  ,  IBUF(40  ) 
DIMENSION  X(ii)O)  ,Y(i00>  ,  IPR  AM  (  S ) ,  IGCB  ( 1 92 )  ,  IOBUF<20)  ,  JCHAR<6) 
DIMENSION  D£LX( 100) ,D£LY< 100) 

DATA  JCHAR/iHO , IH+ , iH* , iHX , iH# , IH$/ 

CALL  RMPAR ( IPRAM ) 

LU»IPRAM<i ) 

LUG=IPRAM< 2) 

IF < LUG . NE . 0 ) GOT 0  IS 
WR ITE(LU, 10 ) 

10  FORMAT (  “  OUTPUT  ON  GRAPHICS  TERMINAL  TYPE  0",/, 
l  "  ON  PLOTTER  TYPE  i " ) 

READ  (  LU,  * ) I LUG 

C  DEFINE  LU  AND  ID  NUMBERS  OF  PLOTTER  AND  GRAPHICS  TERMINAL 

LUG=24 

IF ( ILUG . EQ . 1 )LUG=20 
IS  KK  =  0 

1FILU.LE. 0>LU-i 
ID»i 

IF <  LUG , EQ . 20 ) ID=2 

C  TO  CENTER  CHARACTER  SET  HALF  WIDTH  &  HALF  HEIGHT 

C  DEPENDING  ON  IF  USING  PLOTTER  OR  CRT 

HW=U  .5 
HH=0  .  S 

IF  < ID . EQ . i )  GOTO  20 
HW=Q .333 
HH=0 .25 
20  KK=KI<  +  1 

25  WR ITE ( LU , 30 ) 

30  FORMAT <  "  #**  SOURCE  OF  DATA  ***",/ 

1  "  DISC  FILE  TYPE  0",/ 

2  "  XEYBQADRD  OR  TAPE  TYPE  i - ) 

READ  (LU,*)  IFMT 
IF(IFMT.EQ.O)  GOTO  60 

3S  UR ITE ( LU , 40 ) 

40  FORMAT ( "  ENTER  LU  NUMBER  OF  CTU  OR  KEYBOARD") 

R£A0(LU,*)LUT 
IFCLUT.GT.30)  GOTO  35 
45  WR I TE ( LU , 50 ) 

SO  FORMAT ( "  ENTER  NUMBER  OF  DATA  POINTS  (MAX=100)!,> 

READ(LU,*)NP 

IF ( NP . GT . i 0 0 )  GOTO  45 

DOSS  1=1 , NP 

S3  READ  (LUT,*)  X < 1 ) , Y ( I ) > DELX ( I ) , DELY ( I ) 

GOTO  lOQ 

o0  WRITE(LU,65> 

sS  FORMAT!  "  ENTER  NAME  OF  DATA  FILE") 

READ (LU, 70) NAME 
70  FORMAT (3A2) 
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CALL  OPENUDCB,  IERR, NAME, 3) 

IF (IERR .  LT . 0 ) WR ITE  <  LU , 95) NAME 
73  FOR HAT < “  ERROR  CODE  =  ",I4) 

IFUERR.LT.O)  GO  fO  60 

K*i 

30  DO  82  J=i ,40 

32  IBUF(J)=0 

CALi,  READF ( IDCB , IERR  , IBUF ,  40  , LEN ) 

IF (LEN. EG. -1)  GOTO  85 
IFUERR.LT.  0) GOTO  90 
CALL  CODE 

READ ( IBUF ,  * ) X ( K  > , Y< K ) , DELX < K  > , DELY < K > 

K»K  +  i 
GO  TO  80 
BS  NP*K-i 

CALL  CLOSE  ( IDCB  , IERR > 

IF  < IERR . LT . 0 ) WR ITE (LU, 75 ) IERR 
GOTO  100 

90  WRITE(LU,9S)NAHt 

95  FORMAT ( "  FILE  “ ,3A2,"  DOESN'T  EXIST  OR  IS  ALREADY  OPEN”) 

GOTO  25 

C  IF  DOING  MULTIPLE  PLOTS >  SCALE  DA  I  A  TO  FIRST  PLOT 

100  IF(KX.EQ.l)  GOTO  135 
DOl 30  K  =  1 , NP 
IF< I EXP  >105,115,110 
105  X\K)=X(K>/<iO.*#IABS<IEXP>> 

GOTO  115 

.‘,10  X  ( K  ) *X  ( K  )  * ( 1 0  .  ## IAB3  ( 1EXP  )  ) 

US  IF  (JEXP  >120,130, 125 

120  Y(K)*Y(K>/(10.»*IABS<JEXF>) 

GOTO  130 

1 25  Y ( K  > »Y < K ) # ( i 0 . #*I ABS ( JEXP ) ) 

130  CONTINUE 
i3S  WR IT  E ( LU , 140 ) 

140  p OR MAT ( “  *#*  TYPE  OF  PLOT  #**",/ 

1  "  LINE  PLOT  TYPE  i”,/ 

2  "  SYMBOLS  CONNECTED  WITH  LINES  TYPE  2",/ 

3  ”  SYMBOLS  PLOT  TYPE  3“,) 

READ<i_U,*>  JJ 

r 

C  READ  LABELS  AND  LEGEND 

IF(KK.GT.i)  GOTO  210 
WR I TE( LU, 145) 

145  PORMATC  ENTER  X-LABEL, Y-LABEL  &  TITLE  ON  3  SEPARATE  LINES 
READ ( LU , I5u )  IAX,IAY,IHED 
i.  50  FORMAT  ( 15A2/ ,  1 5A2/ ,  20A2 ) 

C 

C  INITIALIZE  PLOT  PACKAGE  4  DEFINE  VIEWPORT 

C 

CALL  P LOTR ( IGC3 , 1 D , 4 , LUG , IOBUF , 20 ) 

CALL  SETAR ( IGCB,2 . 0 ) 

CALu  V1EWP( IGCB,  0 . ,135.  ,0 .  ,100  .  ) 

n 

w 

C  SCALE  DATA  4  DEFINE  WINDOW 
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n  n  o 


WR I TE  <  LU ,  155 ) 

155  f-  OR  MAT  (  “  **#  ORIGIN  LOCATION  **#",/ 

1  '•  TO  FORCE  ORIGIN  TO  START  AT  (0,0)  TYPE  0",/ 

2  "  TO  LET  SCALE  DETERMINE  ORIGIN  TYPE  i"> 

READ(LU,*)IZERQ 

WR ITE  <.  LU ,  160  ) 

loO  FORMAT ( "  #**  SCALING  X**",/ 

1  “  AUTOMATIC  TYPE  (J "  ,  / 

2  "  MANUAL  TYPE  i " ) 

READ ( LU , * ) ISCAL 

IF ( ISCAL . EQ . 0  ) GOTO  170 
WR I TE ( LU , 165 ) 

165  FORMAT (  "  ENTER  XMIN, XMAX , YMIN , YMAX  "  ) 

READ  ( LU , # ) XMIN , XMAX , YMIN ,  YMAX 
170  CALL  PEN < IGCB , 1 ) 

CALL  SCALE  ( X ,NP , XMIN, XMAX , XTIC , IZERO, IEXP , LU, ISCAL ) 
CALL  SCALE  ( Y , NP , YMIN , YMAX , YTIC , IZERO , JEXP , LU , ISCAL ) 
CALL  WINDW < IGCB ,  0 . , ISO .  ,  0  .  , 1 0 Q  .  > 

C 

C  SPECIFY  CHARACTER  SIZE 

C 

CALL  CSIZE ( IGCB , 3 . ) 

DRAW  AXES  AND  LABELS 

CALL  FXD(IGCB,0) 

SXTIC=ABS  <  XMAX-XMIN ) /XTIC 
SYTIC=ABS ( YMAX-YMIN ) /YTIC 
CALL  MOVE < IGCB, 35. ,1 . ) 

IF < IEXP  NE.0)  GOTO  180 
CALL  LABEL (IGCB) 

WRITE ( LUG , 175 ) IAX 
175  FORMAT (15A2) 

GOTO  190 

180  CALL  LABEL ( IGCB ) 

WRITE (LUG, 185) I AX, IEXP 
185  FORMAT ( 15A2 , "  X  10**" ,12) 

190  CALL  M0VE(IGCB,3. ,12. ) 

CALL  LDIR(IGCB,+1 .57) 

IF< JEXP .NE.0)  GOTO  195 
CALL  LABEL (IGCB) 

WRITE ( LUG , 175  > I AY 
GOTO  200 

195  CALL  LABEL (IGCB) 

WR ITE (LUG, 135) IAY, JEXP 
200  CALL  MOVE (IGCB, 40 . ,90 . ) 

CALL  LDIR ( IGCB , 0 . ) 

CALL  LABEL (IGCB) 

WRITE (LUG, 205) I BED 
205  FORMAT (20A2) 

CALL  VIEWP  <  IGCB,  17  .  ,  120  ,  ,  .1 0  .  ,  30  .  ) 

CaLl.  WINDW ( IGCB, XMIN, XMAX,  YMIN,  YMAX) 

CALl.  LAXES  ( I GCB ,  -SXT I C ,  3 YTIC ,  XM IN ,  YM I N ) 

C 
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•£*  04  ro 


c 

c 

210 


PLOT  DATA  POINTS 


MM=0 

IF< JJ.EQ.2)MM=i 
CALL  LINE ( IGCB, MM ) 

0022 0  K-i , NP 

IFfK.EQ.i.OR. JJ.EQ.3)  GOTO  215 
CALL  DRAW < IOCS, X(K ) , Y<K  ) ) 

C  MOVE  TO  <X,Y) , CENTER  CHAR., PLOT  CHAR., MOVE  "CP”  BACK  TO  <X,Y> 

215  CALL  MOVE < IGCB , X < K ), Y < K  )  ) 

IFCJJ.EQ. i )  GOTO  220 
CALL  CP  LOT  ( IGCB,~HW,~HH,-2> 

CALL  LABEL (IGCB) 

WR I TE ( LUG , 225 ) JCHAR ( KK  > 

CALL  MOVE (IGCB , X( K ) , Y (K  > ) 

220  CONTINUE 
225  FORMAT ( lAi ) 

CALL  PENUP(IGCB) 

CALL  LINE ( IGCB ,  0 ) 

WR  I TE ( LU ,  230  ) 

230  FORMAT ( "  LEAST  SQUARES  LINE  ? " , /,22X , “NO  TYPE  0" ,/, 

1  22X , "YES  TYPE  1") 

READ  ( LU , #  >  ZLINE 
WR I TE ( LU , 235 ) 


233  FORMAT <"  ERROR 

BARS  ?",/, 

i 

NONE 

TYPE 

0", 

2 

X 

TYPE 

i"> 

3 

Y 

TYPE 

2", 

4 

X  &  Y 

TYPE 

3") 

READ <  LU ,♦ ) IBAR 

IF ( ILINE . EQ . 0 )  GOTO  270 


C 

C  CALCULATE  START  AND  END  POINTS  OF  LEAST  SQUARES  LINE 

C 

CALL  LSR EG( X , Y , i , NP , SLOPE , YINT , Si , 32 > 

XST =XMIN 
YST-YINT 

YST®  ( SLOPE*XMIN ) +YINT 
240  IF (YINT . LE . YMAX )  GO  TO  245 

XST=<YMAX-YINT> /SLOPE 
YST* YMAX 

245  IF (YINT . GE . YMIN)  GO  TO  250 
XST* <YMIN-YINT> /SLOPE 
YST® YMIN 
250  XEND®XMAX 

YEND=(SLOPE*XMAX)+YINT 
IF<  YEND . L£ . YMAX  >  GO  TO  255 
XEND»( YMAX -YINT ) /SLOPE 
YEND=YMAX 

255  IF ( YEND. GE. YMIN)  GO  TO  260 
YEND=YMIN 

XEND-(YMIN-YINT) /SLOPE 
C 

C  PLOT  LEAST  SQUARES  LINE 
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260 


IFCILINE.EQ.Q)  GOTO  270 
CALL  MOVEC IGCB,XST, YST) 

CALL  DRAW ( IGCB , XEND , YEND) 

CALL  PENUP  < IGCB ) 

LUX=LU 

XF ( ID . EQ . 2 )  LUX=6 
WR I TEC  LUX , 265  > SLOPE, Si , YINT , S2 
265  FORMAT  (  "  SLOPE* "G12  .  5"  ,  +-,,Gi2.6/"  Y-~j.nT=“Gi2 . 5"  , 
C 

C  COMPUTE  AND  PLOT  ERROR  BARS 

C 

270  IF ( I BAR , EQ , 0 )  GOTO  335 

C  SCALE  DEVIATIONS  BY  VALUE  OF  EXPONENT 

1)030  0  K  =  i  ,  NP 
IFCIEXP >275,285,280 

275  DELXCK)=DELXCK>/CiO.##IABSCIEXP)) 

GOTO  285 

280  DELXcK)=DELX(K)*CiG . **IABS C IEXP )  ) 

285  IF ( JEXP ) 290 ,300, 295 

290  BELY CK)=DELYCK)/CiQ.**IABS< JEXP) ) 

GOTO  300 

295  DELYCK ) =DELY < K ) *C 1 0 . **IABSC JEXP > ) 

300  CONTINUE 
NOBAR=0 

BARX=AB3CXMAX~XMIN)/i 00 . 

BARY-ABS<YMAX-YMIN>/67. 

00325  K*i , NP 
XBAR i =X ( K ) -DELX < K ) 

XBAR2=X(K)+DELXCK> 

YBARi*Y(K) -DELYCK) 

YBAR2=Y ( K ) +DELY  C  K ) 

GOTO ( 3i0 ,320,3i0)  IBAR 
305  GOTO  335 

C  DRAW  X  ERROR  BARS 

C  IF  ERROR  IS  SMALL  DON'T  DRAW  ERROR  BARS 

310  IFCDELX(K) .LE.BARX)  NOBAR=l 
IF(DELXCK) .LE.BARX)  GOTO  3i5 
CALL  MOVE ( I GCB ,XCK),YCK>) 

CALL  CP LOT  c IGCB ,0,3,0, , -2) 

CALL  DRAWC IGCB , XBAR2 , Y  <K  )  ) 

CALL  MGVEC IGCB , XBAR2, Y C K ) +BARY ) 

CALL  DRAWC IGCB, XBAR2, Y CK )-BARY > 

CALL  MOVEC IGCB , X  C  K ) , Y  C  K ) ) 

CALL  CPLQTC IGCB, -0 . 8,0 . ,  -2) 

CALL  DRAWC IGCB , XBAR 1 , Y C K ) ) 

CALL  MOVE  C IGCB , XBAR i , Y  C  K ) +BARY ) 

CALL  DRAWC IGCB ,XBARi , YCK ) -BARY) 

315  IF  (IBAR,  EGJ.  i)  GOTO  32S 
C  DRAW  Y  ERROR  BARS 

320  IFCDELYCK) .LE.BARY)  NOBAR =1 
IF  ( DEi_Y  C  l<  )  .  LE  ,  BAR  Y )  GO")  0  525 
CALL  MOVEC IGCB ,XCK),YCK)) 

CAL.-  CPLOTCIGCB,U . ,0 .5,-2) 

CALL.  DR  AW  C  IGCB  ,  X  (  K  )  ,  YBAR2 ) 

CALL  MOVEC  IGCB ,  X  C  i<  )  -BARX  ,  YBAK2 ) 
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S* 


33  0 
333 

34  U 


C 

343 

350 


360 


CALL  DR AW  <1 GCB , X  <  K ) +BAKX , YBAR  2 ) 

CALL  MODE < IGCB ,X(K),Y(K)) 

CALL  CPLOT(IGCH,0  .  ,-i)  .S,-2> 

CALL  DRAW< IGCB,X(K> ,YBARi) 

CALL  MOVE ( I GCB , X ( K ) -BAR X , YBAR i ) 

CALL  DR  AW ( I GCB , X ( K ) +BAR X , YBAR i ) 

CONTINUE 

CALL  P'ENUP ( I GCB ) 

IF < NOBAR .EQ. i)  WRITE <LU , 330 ) 

FORMAT ( "  *#*  SOME  ERROR  BARS  WERE  TOO  SMALL  TO  PLOT  **#*■) 

WR ITE ( LU  >  340 ) 

FORMAT (  ”  TO  EXIT  TYPE  ()">/> 

X  “  IF  YOU  WANT  ANOTHER  PLOT  ON  SAME  AXIS",/, 

2  "  AND  ALL  X  4  Y  VALUES  ARE  WITHIN  THE",/, 

3  "  SCALES  UF  THE  FIRST  PLOT  TYPE  i "  > 

R£AD(L'J,#> I DONE 

IF ( KK , EU . 6 >  GOTO  360 
SELECT  PEN  COLOR  IF  USING  PLOTTER 
IF  < ID . NE . 2 . OR . IDONE . NE . i )  GOTO  355 
WR ITE ( LU , 350 ) 


FORMAT <  " 

*##  PEN 

COLOR  SELECTION 

II 

BLACK 

TYPE 

i", 

II 

RED 

TYPE 

*  ,  II 
*-  > 

II 

GREEN 

TYPE 

II 

II 

BLLE 

TYPE 

4") 

READ <  LU , # )  IPEN 

IF  (IPEN.LT.i  .OR  .  II 

’EN.GT.4)  GOTO  3 

4S 

CALL  PEN( I GCB, I PEN) 

IF < IDONE. EQ, i)  GOTO  20 
CALL  PEN ( IGCB , 0 ) 

CALL  PLOTR  < IGCB, ID, 0 ) 
END 
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F7N4  ,  L 
C 

ij 

L, 

c 

c 

c 

c 

c 

c 

c 


10 

xS 

c 

20 


2  S 


JO 


40 


43 


3  0 


SUBROUTINE  SCALE  (X ,  NP  ,  XMIN ,  XflAX ,  i  JC,IZERO,  iilXi’ ,  LU ,  ISC 

THIS  ROUTINE  COMPUTES  MAX  4  MIN  VALUES > SCALES  OAT A  I'U 

E  FORMAT  AND  DETERMINES  THE  NUMBER  OF  TICK  MARKS  PER  A 

X  -  ARRAY  1U  BE  SCALED 

NP  -  NUMBER  01-  POINTS  IN  ARRAY  X 

XM  IN  -  MINI NUN  VALUE  OF  X 

XMAX  -  MAXIMUM  VALUE  OF  X 

TIC  -  NUMBER  OF  TICK  MARKS  UN  AXIS 

1ZERO  -  SET  TO  0  TO  FORCE  ORIGIN  TO  ( 0 , 0 ), NORMALLY  =  i 
I EXP  -  EXPONENT  UF  BASE  10  TO  WHICH  X  IS  RAISED 
LU  -  LOGICAL  UNIT  NUMBER  OF  TERMINAL 

I3CAL  -  SCALING,  l)  FOR  AUTOMATIC,  i  FOR  MANUAL 
DIMENSION  XCiQG) 

PIND  MAX  AND  MIN  VALUES  OF  X 
IFCISCAL.EQ.  DGOTO  15 
XMAX=X< 1 ) 

XHIN*X< i ) 

DOS  I =2 , NP 

IF(Xa)  .  GT  .  XMAX  )  XMAX=X  <  I  > 

I F ( X ( I ) . LT . XMI N ) XM IN=X ( I > 

CONTINUE 

1F<I2ER0.NE. O.OR.XMIN.3E. 0)  GUI  0  IS 
WRITE( LU , 10 ) 

FORMAT  < "  **  DATA  HAS  NEGATIVE  VALUES  ORIGIN  CAN  * T  BE  ( 

IZERQ-l 

IEXP=0 

SCALE  DATA  USING  E  FORMAT  AND  SAVING  EXPONENT  (IEXP) 

IF  <  ABS ( XMAX ) , GE . i . OR . ABS (XMIN)  . GE . 1 )  GOTO  30 
DOSS  K=i,NP 
X(K)=X<  10*10  . 

XHIN=XMIN*1 0 . 

XMAX=XMAX#10 . 

IEXP=IEXP+i 
GOTO  20 

IF<A63<XMAX) . LE . lOOO. AND. ABS(XMIN) . LE . 1000)  GOTO  40 

D03S  K -i , NP 

X(K)-X(K)/10. 

XMIN-XilIN/10. 

XMAX=XMAX/10 . 

IEXP-IEXP -1 
GOTO  30 

DETERMINE  INTERVAL  FACTOR 
DIF “ABS ( XMAX-XM1N ) 

IF (DIF, GT.S. )GOTO  4S 
FACTR=1 . 0 
GOTO  80 

IF (DIF , GT . 10 . ) GOTO  50 
FACTR*2. 0 
GOTO  80 

IF ( DIP , GT . 25 . ) GOT 0  SS 
FACTR*S .  IJ 
GOTO  80 

IF ( DIF . GT . SO > GOTO  60 
P'ACTR*ii) .  0 
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GOTO  b'O 

Oii  IF <  DIF  . GT . 125 )  GOTO  63 
FACINGS. 

GOTO  80 

83  IF  ( .OIF  .  GT  .  250  )  GOTO  70 

FACTR=50 . 

GOTO  60 

70  IF <  DIF . GT . 800 ) GOTO  75 

FACfR=iU0 . 

GOTO  30 

75  FAC TR =2 00 . 

80  IF  <  XMIN . GE . 0 , ) GOTO  85 
C  **##  FOR  NEGATIVE  NUMBERS  **** 

MIN=<  XMIN/FACTR  ) -0 . 799999 
XMIN=MIN#FACTR 
IF  ( XMAX.  GE.  0)  GOTO  90 
MAX= ( XiMAX/FACTR ) 

XrtAX=rtAX#FACTR 
GOTO  95 

C  ****  FOR  POSITIVE  NUMBERS  #*#* 

85  MiN=XMIN/FACTR 

XMIN=MIN*FACTR 

9  J  rtAX=  (,  XMAX/FAG TR  )  *  0 . 999999 

XMAX-rtAX#FACTR 

C  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
98  If < IZERQ.EG. 0>  XMIN=0. 

DIF-XMAX-XMIN 

8  DETERMINE  NUMBER  OF  TICK  MARKS  PER  AXIS 

IF  <  DIF . NE . 1 . )GOTQ  100 
110=1 . 

GOTO  130 

100  IF ( DIF . GT , 1 0  >GOTQ  105 

DIF=DIF*1Q , 

GOT O  i 00 

105  IF  <  DIF . LE . 100) GOTO  110 

DIF=DIF/10. 

GOTO  105 
HC  1  IC=3 . 

IF ( DIF . EG . 1 0  0 . >  GOTO  125 
IFCDIi-  .EG. 80.  > GOTO  115 
IF  <  DIF . GE , 60 . ) GO iO  130 
I F  <  D I F . EQ , 5  0 . ) GO  TO  125 
IFOiF.Efl.40.  >COTO  115 
IF t  DIF . EQ , 30 , ) GOTO  130 
IF  <  DIF . EG . 25 . >  GOTO  125 
IF  CD  IF.  EG}.  20  .)  GOTO  120 
IF (DIF. EG. 15. ) GOTO  130 
IF (DIF . EQ . 12 . 5) GOTO  125 
115  TIC=4 . 

GOTO  130 
120  TIC- 2. 

GOTO  130 
i83  i  IC=3 . 

1 50  CONTINUE 
RETURN 

END 
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FTN4 ,  L 

SUBROUTINE  LSREG  < X, Y, IFRM , ITO, SLOPE, YINT , SDSLOP , SDYINT > 
PROGRAM  CALCULATES  LEAST  SQUARES  REGRESSION 
X  -  X  ARRAY 

Y  -  Y  ARRAY 

IFRM  -  STARTING  POINT  OF  INTERVAL 

ITO  -  ENDING  POINT  OF  INTERVAL 

SLOPE  -  RISE  OVER  RUN  OF  L.S.  LINE 

TINT  -  Y  INTERCEPT  OF  L.S.  LINE 

3D3L0P  -  STANDARD  DEVIATION  OF  SLOPE 
SDYINT  -  STANDARD  DEVIATION  OF  Y  INTERCEPT 
DIMENSION  X(100),Y<100) 

FN=ITO-IFRM+i 
TX*  0 
ZY=0 
XY  =  0 
UY=0 
SY  =  0 

DO  100  K=IFRM,ITQ 
SY=SYKY<K)*#2> 

WY-WY+ ( X ( K ) #Y ( K ) > 

XY=XY+X(K) 

ZY=ZY+Y(K) 

100  TX=TX+(X(IO##2) 

TY=<ZY*#2)/FN 
XZY-i XY#ZY )/FN 
XYN=(WY-XZY>**2 
XYD*TX-< <XY**2>/FN> 

SX=<XY#*2>/FN 

SDSLOP=SQRT<  < SY-TY- < XYN/XYD > )/(  (FN-2.  >#<TX-3'X> )  > 
SDYINf=SQRT<  <SDSL0P*#2)*<  TX/FN) > 

5LOFE=<  <FN*WY)-(XY*ZY) )/<  <FN*TX)-(XY**2) ) 
YINT=<ZY-<SLOPC*XY) )/FN 
RETURN 
END 
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